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• ^H , 1. Introduction 

X 

\—{ ' The solar neutrino data [1, 2, 3, 4, 5] and the atmospheric neutrino 

experiments [6, 7, 8, 9, 10] provide strong evidence for neutrino oscilla- 
tions. In the framework of the two flavor neutrinos, these experimental data 
are explained by two sets of the oscillation parameters (A?7ZQ,sin^ 2^0) ~ 
(O(10-5eV2), 0(10-2)) (small angle MSW solution), (©(lO'^eV^), 0(1)) 
(large angle MSW solution), or {O{10~^^eV'^),O{l)) (vacuum oscillation 
solution), and (Am^^^^, sin^ 26iatm) - (10-2•^eV^ 1.0). 

Without loss of generality we assume that |A?tt-2i| < |Am32| < IAtoIj^I 
where Am?- = mf — m?. If both the solar neutrino deficit and the atmo- 
spheric neutrino anomaly are to be solved by energy dependent solutions, 
we have to have Am-li — Awq and Am32 — Am^| , i.e., we have mass 
hierarchy in this case. Therefore I will assume mass hierarchy in the three 
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Fig. 1. Triangle of unit height for 1^3. The vertical position of the state in the 
triangle graph is |C/e3p, which indicates deviation from the two flavor mixings. 

flavor framework throughout this talk. I will adopt the triangle represen- 
tation which has been introduced by Fogli, Lisi, and Scioscia [11]. Fig. 1, 
which will play a role in the analysis of atmospheric neutrino data, rep- 
resents how the most massive state 1^3 mixes with three flavor eigenstates 
with the coefficients |C/a3p {a = e,fi,T), where |f7a3p are the elements of 
the MNS mixing matrix U [12]: 

/ Uel Ue2 Ue3 \ 
U = [ U^l U^2 U^3 ■ (1) 

V f^rl Ur2 UrS J 

2. Constraints from atmospheric neutrino anomaly 

The atmospheric neutrino data of Superkamiokande have been analyzed 
by [13, 14] in the three flavor framework, where smaller mass squared differ- 
ence ATO21 is ignored. The two flavor analysis of the most up-to-date data 
by the Superkamiokande group shows that the allowed region of the mass 
squared difference is 1 x IQ-^eV^ < Am^ < 7 x IQ-^eV^ at 90% CL [8]. The 
analyses in [13, 14] are strictly speaking different from the original one in 
[8], since the full data which are binned with respect to the energy as well 
as the zenith angle are not used in [13, 14]. The analysis in [13] has been 
updated with the recent data for 850 days, where the upward going /j, data 
[9] have also been incorporated. It was found that the region of the mass 
squared difference which is as small as 5 x 10~^eV^ is allowed at 90%CL (cf. 
Fig. 2). 
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Fig. 2. The allowed regions for various Am'^2 by the constraints of atmospheric 
neutrino data of the Superkamiokande contained and upward going /i events, and 
the CHOOZ reactor data. 



On the other hand, the CHOOZ group has updated their result on 
P{i>e -^ i^e) in the reactor disappearance experiment [15], and the mass 
squared difference is limited to Am^ < 7 x 10~*^eV^ for the maximum mix- 
ing. In our three flavor scheme with mass hierarchy, the disappearance 
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Fig. 3. The allowed regions for various Am'^2 by the constraints of atmospheric 
neutrino data of the Kamiokande contained events, the Superkamiokande contained 
and upward going fi events, and the CHOOZ reactor data. All the shadowed regions 
are located near the i^^ — v^ line. 



probability for the CHOOZ experiment is given by 

P{ue -^Ue) = l- sin2 2^i3 sin^ I ^^g^ 



(2) 



so if Am'^2 > 7 X 10 ^eV^ then sin^ 2^13 has to be small. The allowed region 



which is obtained by the constraints of the Superkaniiokande atmospheric 
neutrino data and the CHOOZ data is given in Fig. 2. As is seen in Fig. 2, 
relatively large 6*13 is still allowed for Am'^2 < 1 x 10~^eV^. 

To obtain more stringent bound on sin^ 2^13, I include the three flavor 
analysis [f 6] of the contained events of the Kamiokande atmospheric neu- 
trino data [6], for which the value of the mass squared difference in the 
allowed region tends to be higher than that of Superkamiokande. In fact 
Ar?7.32 < 2 X 10~^eV^ is excluded at 90%CL by including the Kamiokande 
data. Combining the atmospheric neutrino data of Superkamiokande, Kamiokande 
and the CHOOZ reactor data, I have obtained the allowed region which is 
depicted in Fig. 3. Fig. 3 shows that sin^ 2^13 < 0.1 has to be satisfied, 
which is basically the consequence from the CHOOZ data. 

3. A possible way to measure O13 

As we have seen in sect. 2, the data of atmospheric neutrinos and 
the CHOOZ experiment gives sin 2^13 <0.1. The two sets of parameters 

(Ar?^2l,sin^ 2^12) and (A?TT-32,sin^ 2^23) will be determined with more and 
more accuracy in the future by various experiments of solar and atmospheric 
neutrinos, respectively, so the next thing we would like to pursue is to 
determine ^13. There have been discussions on the future intense muon 
beam [17] which could be hundreds times as high as the present one, and it 
would enable us to have very long baseline experiments, where the neutrino 
path length is comparable to the radius of the Earth. In this talk I would 
like to point out that the oscillation probability -P(f^ -^ fe) (in the case 
of Arrigg > 0) or P{T'^ — > J^e) ij^ the case of A?tt-32 < 0) is enhanced for a 
certain region of the neutrino energy due to the matter effect of the Earth 
when sin^ 2^13 >0.01 and therefore it is possible to deduce the magnitude of 
^13 by measuring experimentally P{v^ -^ Ve) or P{i'^ -^ Pg) as a function of 
the neutrino energy in very long baseline experiments which may be possible 
with the intense muon beam technology in the future. 

Let us now consider the situation where A£^2i is completely negligible. 
In that case the positive energy part of the Dirac equation for three flavors 
of neutrinos in matter is given by 

= m\v^\ (3) 

with 

M = [/diag(0,0,A£;32)?7-^ + diag(A,0,0) 




De 



16*23^7 



e^^'i^^Miag (0, 0, AE32) e''"'''-'^" + diag {A, 0, 0) g-^'^^aAT^ 



-J623A7 , 



^g*923A7gi9i^' 'As 



5^diag (1, 0, 1) - ^-diag (1, 0, -1) 



xe-^'^is* 'A5g-ie23A7^-i^ (4) 

where the unit matrix diag(£'2, E2,E2) which contributes only to the overall 
phase has been subtracted from M, AEij = Amfj/2E, E is the neutrino 
energy, A = \/2GFNe{x) stands for the matter effect [18] of the Earth, 
D = diagfe , l,l] and the standard parametrization [19] has been used 
for the MNS matrix (1) 



U 



/ C12C13 S12C13 5136 

-S12C23 - Cl2S23Sl3e*'^ C12C23 - Si2S23Sl3e*'^ S23C13 

\ -512523 - Ci2C23Si3e*'^ -C12S23 - Si2C23Sl3e*'^ C23C13 





= De^^23A7giei3A5^-lgi0l2A2^ (5) 

with the Gell-Mann matrices 

/ -i 
A2 = i 1 , A5 = ( I , A7 = I -i I . (6) 
V 

^13 is the mixing angle in matter given by 

tan2<"'^ Ai^32 sin 2^13 

^^ AE32 cos 2013 ±^ 

as in the two flavor case [18] , and 



B^^^ = V (A^32 cos 2^13 ± Af + (AE32 sin 2^13)^ (8) 

where (+) sign is for antineutrinos, as the sign of A is reversed for an- 
tineutrinos. Assuming the constant density of the matter, the appearance 
probability ^(t'^ -^ Ve) and P{i'^ -^ i>e) in the three flavor framework can 
be written as 

Pi-, - ^e) = sl, sin^ 26^,'^' W (^^) , (9) 

respectively. Note that the only difference between the formulae (9) and 
that in vacuum is sin^ 20j^ which is replaced by sin^ 2^13 in vacuum. 



0.5 
0.4 

> 0.3 

t 
=1 
^ 0.2 

Q. 

0.1 



0) 

> 






0) 

t 
> 



sin 2613 
= 1.2x10"^ 

(613 = ■< °) 



.=£^ 



Itl'-llll 



0.1 



10 




Ev/GeV/(Am32/10"^ev2) 



0.5 
0.4 
0.3 



sin^2ei3 
= 1.1x10"^ 

013 = 3 °) 




Ev/GeV/(Am32/10"^eV^) 



L=1 2800km 
(cos 0=-1 .00) 

L= geookm 

(cos 0=-O.75) 

L= 6400km 
(cos 0=-O.5O) 

L= 3200km 
(cos 0=-O.25) 



Fig. 4. The appearance probability P{v^j, -^ Ve) as a function of Eu//S.m\2 for 
A?7i|2 > 0. A?7i2i — 0, 6'23 = 7r/4 have been assumed. In the case of Avrigj < 
one has to look at P{Vf^ -^ v^) instead. 



Thus the effective mixing angle in matter is enhanced if A?tt-32 > and 
Aii^32 cos 2^13 — A becomes small for some E. On the other hand, if A?7i32 < 



0, then there is no enhancement in the probability Pivn 



but the 



probability -P(^'^ -^ Ve) is enhanced instead. The probabihty ^(t'e -^ ^x) 
with matter effect of the Earth has been discussed by many people in the 
framework of two flavors [20] and three flavors [21]. The equation (9) is 
almost the same as that for two flavor case [20], the only difference being that 
the counterpart to Vf, is the linear combination S23f^ + C23Z^r — ('^^ + ^r)/\/2- 

I have computed numerically the probability P{v^ -^ Ve) for ^23 = 7r/4, 
6)i3 = 1°, 3°, 5°, 7°, 9° (or sin^ 26I13 = 1-2 x lO'^, 1.1 x lO'^, 3.0 x lO'^, 5.9 x 
10-2, 9.5 X 10"2) and for cos 9 = -1, -0.75, -0.5, -0.25 (or L = 12,800 km, 
9,600 km, 6,400 km, 3,200 km) where the zenith angle G, the neutrino path 
length L and the radius R of the Earth are related by L = — 2i?cos 0. The 
results are shown in Fig. 4. The shape of the probability P{i'^ -^ Ve) is 
almost the same as that for the two flavor oscillation [20] , but it is scaled by 
the normalization S23 (cf- (9)). As can be seen in Fig. 4, the case of cos B = 
— 1 (or L = 12,800 km) is most advantageous to get the enhancement in the 
probability for ^13 >2°, while in the case of cos0 = —0.25 it is difficult to 
see the enhancement for smaller values of ^13. 

To measure the probability in practical experiments, one has to measure 
the momentum of the recoiled nucleon as well as that of the outgoing charged 
lepton in a quasi elastic scattering u^ + N -^ i+N'. From Fig. 4 we see that 
the maximum probability is obtained for Ey/GeY^ 1.2 x Am.32/(10-^eV2) 
with COS0 = —1 for each value of ^13. If we assume /S.m^2 — 3.5xlO-^eV2 
which is the best fit value in the Superkamiokande atmospheric neutrino 
data [8], then the maximum probability is obtained for Ey c^ A GeV. The 
smaller | AmggI becomes, the better it works, since the cross section of quasi 
elastic scatterings decreases as the neutrino energy increases [22]. 

4. A possible way to measure 5 

There have been a lot of works which discussed CP violation in neutrino 
oscillations [23, 24]. From the oscillation probability in vacuum 

f /S. EjkL 
\ 2 



P{^a -^ yp; L) = d^p-^Y, Re [UajU^jKkU^k) sin^ ( 

j<k 



+2 J2 Im (UajU^jKkUpk) sin (AE.kL) , (10) 

j<k 



the CP violation in vacuum is given by 

P{l^a -^ 1^13) - P{l>a -^ 1>I3) 

= 4 ^ Im (UajU^jU^f^U^k) sin (AEjkL) 



j<k 
4 J [sin {AE12L) + sin {AE23L) + sin (A£;3iL)] , (11) 



where 

J = lm{UaiU*^,U:,^Up2) (12) 

is the Jarlskog factor, and 

lm[u^,U;,U:2Up2) =lni[u^2U*p2U:,Ups) =lm(u^3U*p^K,Upi) (13) 

has been used. This Jarlskog factor, which is written as 

J = ci3 sin 26*12 sin 26*13 sin 2^23 sin 6 (14) 

in the standard parametrization, contains a small factor sin 2^13 which is 
constrained by the CHOOZ data ( <\/0.1), and a factor sin 2^12 which is 

small in the case of the small mixing angle MSW solution (sin^ 2^12 — 
6 X 10~^). So in general the Jarlskog factor is expected to be small. In 
vacuum the CP violation happens to be the same as the T violation. 

On the other hand, in the presence of matter, the expression (10) for 
the probability is modified. The eigen matrix in matter can be formally 
diagonalized by a unitary matrix V: 

Udi8.g{-AE2i,0,AE32)U-' + dmg{A,0,0) = Vdmg{Ci,C2,C3)V-\{15) 

Assuming constant density, the oscillation probability can be written as 

P{u^ ^ up; L) = 6ap-Aj2 ^^ i^ajV^iKkV^k) sin^ f^^" 

j<k 

+2 J2 Im {VajV^jV:,Vpk) sin (A^.L) , (16) 

j<k 

as in the case of the probability in vacuum. The sign for the matter term A is 
reversed for antineutrinos Ua and the unitary matrix V and the eigenvalues 
Q for Ua are different from V and Cj for neutrinos i^q. Therefore it is not 
illuminating to see the CP violation in matter [24], since it contains terms 
which vanish in the limit 5^0 and those which do not: 



P{ua -^ Vfi) - P{va -^ Vfi) 



-4^ Re 

j<k 

+2^Im 
j<k 



VajV^^V:,Vf3ksm' ( ^^) - VajV^^V:,Vpksm' (^^ 



VajV^,V:,Vpk sin (AQkL) - %jV^,Kk%k sin (AC.feL) , (17) 
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where ^Cjk = Cj ~Ck- It has been pointed out [25] that T violation is useful 
to probe the CP violating phase in the presence of matter. In fact from (16) 
we have T violation in matter: 



AP = P{v^ ^ up) - P{vp ^ z/,) 

= 4^Im(K.,V9',F4y^,)sin(AC,,L) 



j<k 
4 J [sin (ACi2i) + sin (AC23^) + sm (ACsiL)] , (18) 



where 



J ^lui[V^iV^^V:^Vp2) (19) 



is the modified Jarlskog factor in matter and ^Cjk = Cj ~ Ck- Here I will 
consider T violation under two situations where one of the mixing angles in 
(19) is enhanced due to the matter effect of the Earth in our scheme with 
mass hierarchy. 

4.1. (a) \AE2i\ < 1^1 ~ |A^32| 

First let us consider the case where |A£'2i| <C |^| — |A£'32|. This case 
has been discussed by Arafune and Sato [24]. The unitary matrix which 
diagonalizes the eigen matrix (15) to first order in \AE2i/AEs2\ is given by 



V EE e*''23A.g.eff<-'A5exp^iAi5;2i 



— ci3Si2Ci2(Aisin(5 - A2Cos(5) 
1 1 

-Sl3Cl3S?2^5 Sl3Si2Ci2(A6Sin(5 + A7COS(5) 



tt+ — u_ tl+ 



,(20) 



where I have assumed Aii'32 > 0, 9^^ and B^ ' are given in (7) and (8), 
respectively, and 






Ai = 

u± = ^{AEs2 + A±B'^-^), 

Oi3 = Oi3-9fl, , Si3 = sin6'i3, ci3 = cos6ii3. (21) 

The modified Jarlskog factor J'l in this case is given by 

^1 = -^^^^ ^ sin 2012 sin 2^^'"' sin 2^23 sin 5. (22) 

ylcos2C7i3 8 
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If AE32 < then we have to look at AP = P{i'a -^ ^p) — P{^i3 ~^ ^a) 
instead, and 6^ and B^' have to be replaced by 0^ and B^^\ respec- 
tively. Ji is similar to the Jarlskog factor J in vacuum, but one important 
difference between J and Ji is that Ji has sin20j^ which could be en- 
hanced by the matter effect of the Earth. Another difference is that Ji has 
a factor AE21/A whose absolute value is small by assumption and therefore 
the T violating effect under this condition (a) is supposed to be small. 

4.2. (h) |A^2i| ^ \A\ « IAE32I 

Next let us consider the case where |A£'2i| — |A| ^ |A£^32|. In this case 
it is more convenient to adopt the original parametrization for U proposed 
by Kobayashi and Maskawa [26]: 

ci -S1C3 -S1S3 

U = I S1C2 C1C2C3 - 52^36*'^ C1C2S3 + S2C3e*'' 
S1S2 C1S2C3 -I- C2S3e*'' C1S2S3 - C2C3e*'^ 
= e-*^2A7g-i0lA2pgi03A7^ (23) 

where V = diag(l,l,— e j. The eigen matrix (15) is decomposed as the 
zero-th order term Mq and the first order contribution Mi with respect to 

|A^2i/A^32|: 

t/diag (-A^2i, 0, A^32) U-'^ + diag {A, 0, 0) = U{Mq + Mi)t/-^ (24) 



with 



Mo = (0,0,AE32), 

Ml = diag(-A£;2i,0,0)-hC/"Miag(A,0,0)C/ 

= e~^^«^^p-ie*^'i^Miag(t_,t+,0) e~'^'^^Ve'^^^' 

= diag(^i,^2,6) + mAi +r/4A4 + r?6A6, (25) 



where 



1 



t± = 2 ( ^ - ^^21 ± ^iAE2i cos 2^1 - Ay + (AE21 sin 2^1) 
tan2r^ Ai52isin2ei 



^ ^ AE21 cos 26^1 - ^ ' 

7?1 = (t+ - t_)ciSiC3, r/4 = (t+ - t_)ciSiS3, % = (t-S? + t+Ci)c3S3, 

01 = Oi-ef, Si = sin^i, gi=cos^i. (26) 
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Using perturbation theory in | A£'2i /^£'32 1 > it can be shown that the unitary 
matrix 



V = U e*(''-*^5+''6^7)/'^^32g-#f A2 
diagonahzes U{Mq + Mi)U'^ to first order in |AE'2i/A£': 



(27) 



32 1 



U (Mo + Ml) U-'^ = V diag ( ^^^ + 6, ^^^ - 6,6 + A^32 



y^ 



(28) 



where 



h = 



6-6 



tan 2V'i 



M 



2m 



2(t+ - t_)ciSiC3 



6-6 t-cf +t+s2-(t_s2 + t+cf)ci' 



(29) 



In the hmit |A£;2i/A£'32| -^ 0, (27) becomes 

y ~ [/ e-^'^*'^^ 



C^f^el + S^Ue2 —S^Uel + C^C/e2 Ue3 
CtpU^l + S^t/^2 —Si/jU^i + C^C/^2 ^^3 
C^f^rl + S^Ur2 —S^Url + C^Ut-2 Ur3 



(30) 



where s^ = sin -;/'{, c^ = cos'f/'i • Notice that Vi survives even in the 
hmit I A£'2i / ^£^32 1 -^ 0. Therefore the modified Jarlskog factor J2 is given 
by 



J2 



Im [{CjpUel + S^t/e2)(-S^f^el + C^f^e2) 
X (c^f/^1 + S^f/^2)*(-.S^f^Ml + C^U^2)] 



i.M 



sin 2V'r (1 - U^s) - UelUe2 COS 2i)f 



M 



C/e3 sin 26*2 sin (5. (31) 



In the case of the small mixing angle MSW solution for which we may 
neglect small factors U^^, ^^d Ue2, we have 

J2 = J sin 2^f sin 26'2[/e3 sin 5 ~ - sin 29^ sin 292Ue3 sin 5, (32) 

where we have used the fact tan 2'tpf^ ~ tan 26i^ which follows from (29) 
when 1^3! ^ 1. (32) could be large since J2 does not contain the suppression 
factor |A£'2i/Ai?32| like in J'l, sin20i could be enhanced by the matter 
effect of the Earth (cf. (26)), and the only factor which is always small is 
Ue3 (|[/e3| < VO.1/2). Note that we know A£^2i > from the solar neutrino 
deficit. 
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Fig. 5. T violation P{Vf, -^ Ue) - P{i^e 



Vi, 



in matter (solid lines) and in vacuum 



(dashed lines) for the large (Amji — 1.8 x 10 ^eV , sin 29i2 = 0.76) and small 
(Amii = 5.4x IQ-^eV^, sin^ 29i2 = 0.006) angle MSW solutions, sin^ 29^ = 0.036 
and the set of the best fit parameters (Am§2,sin^ 2^23) = (3.5 x 10^'^cV^, 1.0) for 
the atmospheric neutrino data have been taken as reference values, maximum T 
violation S = tt/2 has been assumed and rapid oscillations due to larger eigenvalues 
have been averaged over. 



4.3. Numerical analysis 

The results by numerical calculations are given in Fig. 5 where rapid 
oscillations coming from the larger eigenvalues are averaged over, maximal 
T violation is assumed (5 = tt/2) and the best fit values [27] of the large 
and small angle MSW solutions have been chosen for the set of parameters 
(Ar7T,2i,sin 2^12), respectively. In the case of the large mixing angle MSW 
solution there is some resonance for Eiy ~ several GeV and Ei, ~ 0(0.1) 
GeV which satisfy the conditions (a) and (b), respectively. In the case of 
the small mixing angle MSW solution, the enhancement is hardly visible 
for the energy E^, ~ several GeV which satisfies the condition (a) but it 
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is significant for Eu ~ 0(10) MeV wfiich satisfies the condition (b). We 
see from Fig. 5 tliat T violating effects are significant for relatively low 
neutrino energy, i.e., E^ ~ 0(0.1) GeV and E^ ~ 0(10) MeV in the case of 
the large and small angle MSW solutions, respectively. For the former case, 
the enhancement around E^, ~ several GeV might be observable at neutrino 
factories [28], where intense beams of Ug and u^ (or i>e and P^) are produced. 
For the latter case, it would be very difficult to see these effects unless we 
can produce intense beams of neutrinos with energy Ei, ~ 0(10) MeV. It 
should be also mentioned that it is hopeless to see T violation effects in the 
case of the vacuum oscillation solutions, since Am^i is extremely small. 

5. Conclusions 

To summarize, I have shown that it is possible to probe the small value 
of ^13 down to sin^ 2^13 >0.01 in very long baseline experiments by looking 

for Ve (in the case of Am'^2 > 0) o^ ^e (in the case of Am32 < 0) appearance 
which is enhanced due to the matter effect of the Earth. If we can produce 
very intense beams of i^^ and Vg with low energy then we may be able to 
determine the CP violating phase S by looking at T violating effects in the 
case of the MSW solutions. 
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